Unconventional proximity effect and inverse spin-switch behavior in a model 
manganite-cuprate-manganite trilayer system 



Juan Salafranca^'^ and Satoshi Okamoto^ 

^Department of Physics and Astronomy, The University of Tennessee, Knoxville, Tennessee 37996, USA 
^Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 

(Dated; today) 

The proximity effect in a model manganite-cuprate system is investigated theoretically. We con- 
sider a situation in which spin-polarized electrons in manganite layers antiferromagnetically couple 
with electrons in cuprate layers as observed experimentally. The effect of the interfacial magnetic 
coupling is found to be much stronger than the injection of spin-polarized electrons into the cuprate 
region. As a result, the superconducting transition temperature depends on the thickness of the 
cuprate layer significantly. Since the magnetic coupling creates negative spin polarization, an applied 
magnetic field and the negative polarization compete resulting in the inverse spin-switch behavior 
where the superconducting transition temperature is increased by applying a magnetic field. 
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Transition-metal oxides have been providing intriguing 
phenomena due to strong electron-electron or electron- 
lattice interactions, such as high-Tc superconductivity 
(SC) in cuprates and the colossal magnetoresistance ef- 
fect of manganites [H. Recent developments in fabri- 
cating atomically controlled heterostructures comprised 
of different transition-metal oxides allows us to explore 
further exotic phenomena that are not realized in bulk 
systems d-Q]. Heterostructures involving cuprates and 
manganites have attracted much attention because of the 
competition between SC and nearly full spin polarization, 
and their potential application as spintronic devices 0- 

This growing interest has rendered manganite-cuprate 
heterostructures the paradigmatic example of multilayers 
composed of oxides with competing ordered states. The 
electronic structure near the interface presents remark- 
able differences as compared to the bulk. The electronic 
charge is redistributed [^, [lH, and orbital reconstruc- 
tion takes place As a consequence, spin polarization 
is induced in the cuprate side [10[. The detailed char- 
acterization of the interface has been complemented by 
systematic studies of the collective properties. 

It has been revealed that the SC critical temperature 
Tc of YBa2Cu307_5 (YBCO) is more strongly suppressed 
when combined in superlattices with ferromagnetic (FM) 
Lai_^Ca^Mn03 (LCMO) Q than with a nonmagnetic 
cuprate [l^l- Even 5 unit-cell (u.c.) thick YBCO can 
become nonsuperconducting depending on the thickness 
of the LCMO layers . Since the c-axis coherence length 
in YBCO is about 1 u.c, these experiments indicate the 
existence of an unconventional proximity effect. 

More recently, a surprisingly robust inverse spin-switch 
effect (ISSE) was discovered in LCMO/YBCO/LCMO 
trilayer systems [l3 , 14 1 . In sharp contrast to the conven- 
tional exclusion between SC and magnetism [l5| . SC is 
favored by parallel alignment of the magnetization in the 
FM layers under an applied magnetic field. Changes in Tc 



in these systems were found to be as high as 1.6 K jl4|. 
A similar ISSE has also been reported in conventional 
ferromagnet/superconductor heterostructures. However, 
the changes in Tc were much smaller (~ 10 mK) (l6. 17|. 
and these are most likely due to stray fields [l^ . A large 
band splitting between the majority and minority elec- 
trons could have such an effect [19|. But, the ISSE is 
shown to disappear when the minority band is above the 
Fermi level, and the mechanism of Ref. [l^ cannot ac- 
count for the ISSE in manganite-based heterostructures. 

These experimental findings indicate that the inter- 
facial phenomena in manganite-cuprate systems lie far 
outside of conventional theoretical models which m erely 
consider the transfer of electrons at interfaces [ll|: Il9| . 
Identifying the proper interfacial interactions and pro- 
viding physical pictures is, therefore, desired not only for 
understanding the experimental results but also for their 
potential device applications. 

-manganite trilayers. As the key ingredient, we con- 
sider the antiferromagnetic (AF) coupling at the inter- 
face confirmed experimentally [lO| . The AF coupling in- 
duces the negative spin polarization inside the cuprate 
region and infiuences SC more strongly than injecting 
spin-polarized quasiparticles without the coupling. We 
found that Tc of such systems is drastically suppressed 
when the thickness of the cuprate layer is reduced, in 
accordance with the experimental reports. The bal- 
ance between the effective field due to the interfacial AF 
coupling and external applied fields naturally explains 
the ISSE. The model also reproduces the characteristic 
length scales observed experimentally [HI and the expo- 
nential decay of the ISSE with cuprate-layer thickness 
[isj . Our work quantitatively links the details of the 
coupling at the manganite-cuprate interfaces and super- 
conducting properties of the multilayers, and offers a co- 
herent picture to understand most experimental results 
reported in manganite-cuprate heterostructures. 

Model. — We consider a [001] trilayer system in which 
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FIG. 1. Model, (a) Realistic manganite-cuprate interface 
and (b) simplified interface. Manganite conduction electrons 
and cuprate 3z^ — electrons are coupled via the antiferro- 
magnetic exchange J'af, while cuprate 3z^ — and — y'^ 
electrons are coupled ferromagnetically via the Hund coupling 
J'lj. Since the Hund coupling is larger than J'af^ the mag- 
netic interaction can be approximated as a direct coupling 
Jaf as shown in (b). Trilayers with (c) an antiparallel config- 
uration and (d) [(e)] a parallel configuration (with an applied 
magnetic field H). H mostly affects the tails of the induced 
moment. 



manganites and cuprates are stacked along the z direc- 
tion. The manganite region is described by the double- 
exchange model Hde with single conduction band for 
simplicity. The cuprate region is described by the BCS 
model HbcS: with the pairing involving Cu dx2-y2 or- 
bitals. At a manganite-cuprate interface, Cu c?322_r2 or- 
bitals become electrically active due to the "orbital recon- 
struction" 0,[iS[2l|. Thed3^2_^2 orbital hybridizes with 
manganese Cg orbitals, and the resulting AF coupling be- 
tween them is J'j^p- Further, the d^^^-r'^ orbital and the 
Cu dx^-y^ orbital on the same Cu site are expected to 
be coupled ferromagnetically via the Hund coupling 
[see Fig. [B^a)]. Since » J'af' the Cu d^^-y^ spin 
is expected to be slaved to the Cu d^^^-r^ spin. Such 
a situation is modeled by a direct AF coupling Jaf be- 
tween a manganese spin and a Cu d,j.2_y2 spin as shown 
in Fig. [TJb). The coupling constant Jaf is the same or- 
der of magnitude as J'^^p but somewhat reduced. Thus, 
the Hamiltonian of our model trilayer system is given by 

HfSF = J2a=R.,L{^DE,a + Hint,a} + HbCS with 

Hde,o.=- t (4c,- + h.c}j - Jh s**'" ■ S,. (1) 

+5 Y ['^h4+fi'^^+r\dii + 4A+f^d.i+ridi^), (2) 
H,nt,c.=Y{{'''^^ldoo + /i.e.) + JAFst'" • "}. (3) 

Here, a=R or L indicates the manganite layer on the 
right or left, Jh the Hund coupling in manganite re- 



gions. c{d) is an electron annihilation operator in a man- 
ganite (cuprate) region, and s^^" = \ Y.aa' 4a'^<y<y"^^'y' ^ 
~ ^J2aa' 4a^<^<^'dicr' with T the Pauli matriccs, 
and Si a localized t2g spin in the manganite regions 
(15*^1=1). The transfer intensity in a manganite (cuprate) 
region is given by t{tij). The hybridization strength at a 
manganite-cuprate interface is given by Va with (ij)' in 
Hint,a indicating the summation constrained for nearest- 
neighbor bonds across the interface, g is the pairing cou- 
pling constant, x{y) is the unit vector along the x{y) 
direction, and Se indicates the band mismatch between 
the manganite and cuprate regions. Considering the ex- 
perimental set up in Refs. [l3, 14 1, i.e., the external mag- 
netic field being applied in in-plane directions, the effect 
of the magnetic field H is introduced by the Zeemann 
term — geS^H with ge = 2, and i running in both Cu 
and Mn regions. 

Since we are interested in the superconducting proper- 
ties, we focus on the temperature range far below the FM 
Curie temperature of manganites. t2g spins in manganite 
regions are then treated as classical degrees of freedom 
with parallel (P) or antiparallel (AP) spin alignment be- 
tween the R and L manganite layers. We then apply 
the Bogoliubov-Hartree-Fock approximation by introduc- 
ing the order parameters nf^'^ = {c\^Cic) ,nfj' = (dl^di^), 
and A,^ = (di^di+ooi) = -{di\di+yi) (d-wave symmetry). 
These order parameters are determined by solving the 
self-consistent equations numerically. 

Each system in our simulations is characterized by the 
thickness of the cuprate middle layer iV, with the lattice 
constant taken to be unity. The total number of unit 
cells along z in the manganite layers is made equal to N 
to reduce the parameter space. In-plane hopping in the 
Cu region and isotropic hopping in the Mn regions are of 
the same order of magnitude and are larger than hopping 
along z in the Cu region and across the interface. Jaf 
and g are derived from the superexchange-type processes. 
The effect of JAF{g) should be enhanced (suppressed) by 
the e-e repulsion, although it is not included explicitly. 
Therefore, we use the following parameters: Wfl = vl = 
0.2i, Jaf = 0.6i, Jh = lOi, g = 0.2t, Uj = t (0.2t) for a 
nearest-neighbor transfer along the xy (z) direction, 5e = 
l.2t, and the mean carrier density in the bulk cuprate 
(manganite) 0.8 (0.7). The bulk chemical potential for 
the cuprate region is about l.5t lower than that in the 
manganite region. This results in a small electron doping 
of the cuprate region in the multilayer model systems 
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As YBCO samples are always in the nearly optimal 
or underdoped sides of the phase diagram, these doping 
effects will enhance the suppression of Tc in thin het- 
erostructures, bringing the system closer to the Mott in- 
sulating phase. Enhanced phase fluctuations in thin het- 
erostructures should also reduce Tc- But these effects are 
not included in the present treatment. The penetration 
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FIG. 2. Tc as a function of the Cu layer thickness A'' in 
the P configuration. Jaf = 0.6t reduces Tc by inducing the 
spin polarization in the Cu region, while Jaf ~ shows the 
weak A'' dependence of Tc demonstrating that the polarized 
quasiparticle injection is insignificant. Inset: Magnetization 
m = rif — ni (circles) and charge density n = n^+n^ (squares) 
profiles for the N — 8 trilayer with Jaf = 0.6t. Ua- = n^"'^^"' 
in the cuprate (manganite) region at 5 < z < 12 (z < 5 and 
z > 12). Results at r = 0{Tc) are shown as open (filled) 
symbols. 

depth of the SC order parameter in the FM region is ex- 
tremely small due to the high spin polarization, and the 
c-axis SC coherence length is short in YBCO, therefore, 
results for trilayers in the P configuration with N > 2 
are representative of experiments in multilayers. 

A larger Jaf implies that the AF coupling has a 
larger effect suppressing the superconducting order, and 
a larger reduces Tc further, but our numerical re- 
sults do not depend on the choice of these parameters in 
a significant way. For different system sizes, A is calcu- 
lated as a function of temperature and the magnetic field. 
Tc is then defined as A{Tc) = 10^2^(0). The numerical 
error is smaller than the point size presented below. 

Unconventional proximity effect. — Figure [2] (main 
panel) shows Tc as a function of N for a trilayer in the P 
configuration. Tc is strongly suppressed by reducing A'^, 
and SC is absent for A^ = 2. Both the thickness at which 
SC disappears, and the length scale for a reduction in 
Tc are very similar to the experimental results Q. For 
Jaf = 0, Tc remains relatively unchanged as a function 
of A^. This demonstrates that the AF coupling plays a 
key role in the phenomena discussed here, while the in- 
jection of highly polarized quasiparticles is insignificant. 

The magnetization density m profiles in Fig. [2] (inset) 
show the origin of the characteristic length scale of the 
unconventional proximity effect, m in the cuprate layer 
is more than 1 order of magnitude smaller than in the 
manganite layer, even for the Cu ions closest to the sur- 
face. In the absence of SC (T = Tc), m decays faster 
into the cuprate layer. Screening sets the length scale for 
the unconventional proximity effect. When a strong SC 
phase is present, it reduces m near the surface, but the 
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FIG. 3. Tc as a function of the applied magnetic field H for 
different N. Circles (triangles) indicate Tc's for the P (AP) 
configuration (see Fig. [T|). Continuous lines are quadratic fits; 
the linear coefficients determine the behavior of Tc{H) for the 
experimentally relevant fields. A wide range of magnetic fields 
is considered (note that t « 0.2-0.5 eV). 



screening takes place in an oscillatory manner. 

ISSE. — Next we consider the effect of an applied mag- 
netic field in trilayer systems, i.e., the ISSE. Figure[3]dis- 
plays Tc as a function of H for the various A^ indicated. 
The filled circles are numerical results, and the lines are 
quadratic fits that will be discussed later. Changes in en- 
ergy associated with the i/'s are small, and it is expected 
that the relevant properties can be expanded in power se- 
ries of the field. For the AP configuration, the system is 
symmetric against the inversion of H , and therefore the 
linear term in the expansion of Tc is zero. For the P 
configuration, the linear term is always positive for the 
different N investigated. This explains the surprisingly 
low field negative magnetoresistance experimentally ob- 
served for trilayers in the P configuration [l3| . 

A qualitative picture of these effects can be drawn from 
Figs. [Hand [21 Without iJ, there are some small induced 
magnetic moments in the cuprate region [Figs. [Ijc) and 
[IJd), Fig. [2] (inset)]. For the P configuration, applied 
fields compensate the effective field due to the AF cou- 
pling with manganites, especially at the exponential tail 
part [Fig. [TJe)]. Thus, Tc increases at low and moderate 
field regimes. In the AP configuration, effective and ap- 
plied fields compensate only at one interface region but 
add up at the other. The ISSE is therefore caused by 
the different effects of the external fields in the P and 
AP configurations. Experimental details, such as shape 
anisotropy, determine the switching fields and would in- 
duce rather nonlinear and hysteretic behavior. 

For the A^ = 6 trilayer with no applied field, Tc in the P 
configuration is lower than in the AP configuration. The 
effective negative field is negligible far enough into the 
cuprate layer, but large near the interfaces. For small A^, 
effective fields from the different interfaces overlap. They 
partially cancel in the AP configuration, while they add 
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FIG. 4. Change in the critical temperature ATc = ~ T^^ 
due to the different relative orientations of the magnetization 
of the manganite layers at _ff = 2 x 10~^t as a function of A'' 
(squares). The exponential behavior is due to the linear de- 
pendence of (H). For the fields accessible to experiments, 
this is the only relevant term, and the linear coefficient b also 
decays exponentially with N (circles). 

up in the P configuration. Therefore, Tc is lower for this 
last case. Nevertheless, nonzero H produces a spectacu- 
lar rise in T^, as the interfacial effective field is very large. 
Even for N = 2, becomes finite aX H > 0.035. The 
sensitivity of the thinnest layers to the magnetization 
orientation might contribute to the reported dependence 
of the SC properties on the manganite layer thickness 
0. But, additional effects such as phase fluctuations not 
considered here might be important in this limit. 

For iV = 8 and thicker trilayers, Tc{H = 0) for the 
P and AP configurations are the same within our nu- 
merical precision. The effective fields corresponding to 
different interfaces do not overlap. The P configuration 
has a higher Tc under an applied field due to the par- 
tial cancellation of the external magnetic field and the 
effective field from the interface, giving rise to the ISSE. 

Finally, we discuss the change in Tc due to the differ- 
ent alignment of the FM layers AT^ = T^ - T^^ as a 
function of thickness. T^f and T^^ change with an ap- 
plied field (Fig. In experiments, the switching fields 
used to control the magnetic configuration depend on 
the details of the samples and the experimental setup. 
In Fig. H ATc is defined as T^ - T^^ at H = 2x IQ-H. 
With these switching field values, ATc decays exponen- 
tially with N (squares in FigH]), in agreement with Ref. 
[isj . The choice of switching fields does not alter the 
exponential dependence of ATc if these fields are small. 

For the fields accessible to experiments, T^ increases 
linearly with the field, while T^^ only changes with the 
second power. The zero field value of Tc is equal for both 
configurations except for the thinnest layers. Thus, for 
the relevant switching fields H^, ATc = bHs + 0{Hf), 
where b is the (thickness dependent) linear coefficient in 
Tf{H). b is plotted in Fig.ll as a function of the cuprate 
layer thickness (circles). It is clear that the change in b 
with respect to the cuprate thickness is the origin of the 



exponential decay in ATc, and that the switching field 
only appears as a vertical shift in the logarithmic scale of 
Fig. m The exponential decay length for the parameters 
used in this Letter is 1=2.7 u.c, whic is smaller than 13 
nm (« 10 u.c.) found experimentally [l3|. Fine-tuning of 
the parameters to make the experimental and theoretical 
results overlap is indeed possible, but vortices (not ac- 
counted for in the present mean field treatment) reduce 
the diamagnetic effect of the superconducting phase and 
are expected to increase the characteristic length of the 
ISSE. Our model provides a new explanation for this ef- 
fect, reproduces the correct exponential behavior, and 
gives a qualitative agreement with experiments. 

To summarize, we investigated the proximity effect 
in manganite-cuprate systems with the AF coupling be- 
tween the manganite and cuprate regions. This coupling 
induces negative spin polarization in the cuprate region, 
in contrast to the positive spin polarization expected from 
a simple charge transfer picture. The effect of the AF 
coupling was found to be much stronger than inject- 
ing spin-polarized quasiparticles. Various experimental 
anomalies are semiquantitatively reproduced based on 
our model calculations, such as the strong thickness de- 
pendence of Tc in manganite-cuprate superlattices and 
the inverse spin-switch behavior of manganite-cuprate- 
manganite trilayers. As in this study, identifying proper 
interactions is crucial to understand the novel phenom- 
ena observed in many other oxide heterostructures. 
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